Negative diffraction pattern dynamics in nonlinear cavities with left-handed materials 
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We study a ring cavity filled with a slab of a right-handed material and a slab of a left-handed 
material. Both layers are assumed to be nonlinear Kerr media. First, we derive a model for the 
propagation of light in a left-handed material. By constructing a mean-field model, we show that 
the sign of diffraction can be made either positive or negative in this resonator, depending on 
the thicknesses of the layers. Subsequently, we demonstrate that the dynamical behavior of the 
modulation instability is strongly affected by the sign of the diffraction coefficient. Finally, we 
study the dissipative structures in this resonator and reveal the predominance of a two-dimensional 
up-switching process over the formation of spatially periodic structures, leading to the truncation 
of the homogeneous hysteresis cycle. 

PACS numbers: 41.20. Jb,42.70.Qs 



I. INTRODUCTION 

In recent years, considerable progress has been realized 
in the understanding of dissipative structures in spatially 
extended systems. The concept of dissipative structures 
in systems far from equilibrium has been first introduced 
in the context of reaction diffusion systems [1, 2]. Their 
formation is attributed to the balance between nonlin- 
earities (chemical reaction or light- matter interaction), 
transport processes (diffusion and/or diffraction), and 
dissipation. Driven optical cavities filled with nonlinear 
media - such as Kerr materials, liquid crystals, semi- 
conductors, or quadratic media - belong to this field of 
research and constitute a basic configuration in nonlin- 
ear optics. It has been shown, both experimentally and 
numerically, that these simple devices support station- 
ary dissipative structures that can be either periodic or 
localized in space (for an overview, see [3-7]). 

Along another line of research, theoretical and experi- 
mental studies have shown the possibility of making left- 
handed materials (lhm), i.e., materials with simultane- 
ously negative permittivity and permeability. Although 
the early lhms were build for microwaves, the first opti- 
cal LHMs have been demonstrated recently [8-14]. In his 
seminal work, Veselago has shown the possibility of light 
propagation in these materials and predicted a number 
of exotic effects like negative refraction, inverse Doppler 
shifts and radiation tension [15]. Several experiments 
have confirmed that LHMs must be characterized by a 
negative index of refraction (n < 0) and, hence, by a 
negative wavenumber (fco < 0) [16, 17]. Linear optical 
resonators filled with LHMs were shown to exhibit unique 
properties, most noteworthy their unconventional disper- 
sion relation [18]. Also, periodic heterostructures that 
combine positive and negative refractive indexes present 
analogies with cavities and show interesting behavior. In 
particular, a new bandgap corresponding to zero aver- 
aged refractive index appears [19, 20]. Despite the im- 
portance of nonlinear phenomena to the propagation of 



radiation in lhms, only a few works have addressed their 
nonlinear properties [21-26]. In particular, it was shown 
that LHMs exhibit a bistable behavior with respect to 
power, which suggests to use them in power limiters and 
all-optical switches. 

To the best of our knowledge, the latter studies ad- 
dressing the nonlinear behavior of lhms consider the 
laser beam as homogeneous in the plane transverse to 
the beam propagation axis and therefore neglect diffrac- 
tion. The purpose of this article is to report on analytical 
and numerical investigations of an externally driven ring 
cavity filled with two materials having indices of refrac- 
tion of opposite signs, taking diffraction into account (see 
Fig. 1). In a very recent study [27] a 2-layer (lhm+RHM) 
nonlinear Fabry-Perot interferometer is considered, but 
the propagation in the lhm is assumed to be linear and 
the study is focused on the behavior of the device when 
the input angle is varied. In that case, surface waves at 
the interfaces are to be considered, which is not necessary 
in our study, where the pump beam is orthogonal to the 
entrance surface. 
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FIG. 1: Geometry of the ring cavity filled with right-handed 
(rhm) and left-handed (lhm) materials. The cavity is driven 
by the external field Ei„. The internal forward and backward 
fields are denoted by Ef and E^, respectively. 

We show that the dynamical behavior of the modu- 
lational instability (mi) is strongly affected by the left- 
handed material layer. This is due to the possibility 
of controlling the magnitude and even the sign of the 
diffraction coefficient by varying the relative thicknesses 
of both layers in the resonator. Finally, we describe 
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the two-dimensional up-switching process that dominates 
the spatiotcmporal behavior of the nonlinear resonator. 
Note that the term negative diffraction has been intro- 
duced in other systems too, such as atomic Bose-Einstein 
condensates [28] and discrete systems [29, 30]. Further- 
more, the possibility of diffraction engineering using pho- 
tonic crystals has been raised in [31], where it was shown 
that the width of solitons may become as small as a few 
times the wavelength of the modulation period. We want 
to propose a completely different mechanism for tuning 
the diffraction coefficient. In addition, the solitons of 
Ref. [31] were studied with only one transverse dimen- 
sion, while here, we consider 1-D and 2-D systems and we 
show that beyond the similarity of their evolution equa- 
tions, the one- and two-dimensional cases exhibit very 
different dynamics. 

In this article, we consider an optical ring cavity, 
which is driven by a coherent light beam of frequency 
v a = cj /27t. The cavity is filled with two adjacent lay- 
ers containing a right-handed material (rhm) and a left- 
handed material (lhm), respectively. We will assume 
that both materials are local, dispersive and have a Kerr 
type nonlinearity in their dielectric behavior. These as- 
sumptions are quite common for traditional (rhm) ma- 
terials, but need some explanation for lhms. First, the 
absence of nonlocality can be justified from the fact that 
lhms are subwavelength-structured materials; coupling 
between resonators over distances longer than the wave- 
length will therefore be negligible with respect to diffrac- 
tion. Second, lhms have been shown to be dispersive 
in both their dielectric and magnetic response [15]. The 
third assumption seems to contradict the magnetic non- 
linearity found by Zharov et. al. [21]. However, their 
nonlinearity is a consequence of a nonlinear dielectric 
inserted between the gaps of the split-ring resonators 
in Smith's metamaterial, forming nonlinear capacitors. 
This will be unlikely the case for optical lhms. If the 
nonlinearity can be introduced by other means, e.g. by 
using nonlinear resonators interacting with the electric 
field [32], the nonlinearity can indeed be modelled by a 
third-order susceptibility for a centrosymmetric medium. 
Furthermore, it can be inferred from the analysis pre- 
sented below that a Kerr nonlinearity in the RHM suffices 
to observe the effects reported in this paper and that the 
lhm may even be simply linear. 

The structure of this paper is as follows. In Sect. II, 
we will derive a nonlinear Schrodinger equation describ- 
ing the propagation of light through the lhm. Subse- 
quently, in Sect. Ill, a mean-field model for the nonlinear 
resonator will be constructed. A linear stability analy- 
sis of this model will be carried out in Sect. IV and the 
dissipative structures emerging due to the modulational 
instability will be discussed in Sect. V. 

II. LIGHT PROPAGATION IN A LHM 

Light coupled into the resonator depicted in Fig. 1 
will propagate through the right-handed layer as well as 



through the left-handed layer. A physical model of this 
resonator should therefore include a description of both 
interactions. 

On the one hand, it is well-known that propagation 
through a right-handed Kerr type nonlinear material 
can be properly treated with the standard nonlinear 
Schrodinger equation (nlse) [33], which gives the enve- 
lope A of the electric field as a function of time t, the 
longitudinal coordinate z, and the transverse coordinates 
x and y: 

where £ is the longitudinal coordinate of a reference sys- 
tem traveling with the group velocity along the cavity, ujq 
is the central pulsation of the pulse envelope, n is the 

(3) 

refractive index of the material at loq, \o its third-order 
nonlinear susceptibility at too, and c is the speed of light 
in vacuum (eoMo c 2 = !)• 

On the other hand, propagation through a nonlinear 
left-handed material has not been considered before. It 
is not obvious at all that the nlse is still valid for left- 
handed materials, because it is commonly derived for 
nonmagnetic materials with [i r — n/no = 1. There- 
fore, in this Section we will outline how the nlse can 
be generalized for dispersive magnetic media. A detailed 
treatment can be found in Ref. [34]. 

In the first step, we model the material behavior by 
the following constitutive equations relating the displace- 
ment field D and the electric field E on the one hand, and 
between the magnetic field H and the induction field B 
on the other hand: 

/+oo 
£ r (i-T)E(r,r)dT + p( 3 \ (2) 
-CO 

" (S, =-///_>"»(*-r„ ( -r 2 ,*-. a) p) 
E(r, n)E{r, r 2 )E(r, T 3 )dT 1 dT 2 dT 3 , 

/+oo 
M r (i-r)H(r,r)dr. (4) 
-oo 

These equations account for the dispersive dielectric and 
magnetic behavior by the convolution integrals and for 
the Kerr nonlinearity by the third-order polarization field 
P^ 3 ). The equations are valid for dispersive magnetic 
materials that are cubic and centrosymmetric. Combin- 
ing Eqs. (2)- (4) with Maxwell's equations, we obtain the 
wave equation governing the evolution of the electric field 
in a nonlinear dispersive magnetic dielectric: 
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^ /■oo ^2 /■oo <-oo 

-V x V x E(r,i) - — / Vrit-r)— e r (r - u) E(r, u) dudr = ,u / /U r (i- 
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r) J^P( 3 >(r,r)dT. (5) 
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One can now proceed with the derivation of an amplitude 
equation for quasi-monochromatic light by invoking the 
slowly-varying envelope approximation and by applying 
the classical multiple scales perturbation technique [33]. 
In this derivation, we represent the coherent light beam 
by the following wavepacket 



E = e4(r,()e i < k »'"» i ' +c.c, 



(6) 



and we obtain a generalized version of Eq. (1) that is 
valid for both right-handed and left-handed materials: 



.8A + c 
d( 2wo7io 
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2 dt 2 2ct/o 
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0, (7) 



with r/o — yje r / Hr, the characteristic impedance of the 
medium. We thus find that the evolution of the envelope 
of the electric field is still described by a nlse, but the 
coefficients of Eq. (7) differ from those of Eq. (1). In 
particular, the refractive index no in the nonlinear term 
of the original equation must be replaced by the charac- 
teristic impedance r] . From our detailed analysis, it is 
now also clear that the coefficient of the diffraction term 
is still inversely proportional to the refractive index, and 
therefore negative in a LHM. 



III. RESONATOR MODEL 



scale t = \/2nt/ {FT), where T is the cavity round- 
trip time and T its finesse. The transverse space co- 
ordinates x and y and the longitudinal space coordi- 
nate z are left unchanged. With these conventions, the 
field circulating in the ring cavity can be described by 
Ef = A(x, y, t) exp [i (koz — coot)], with ko the wavenum- 
bcr. The evolution of the envelope A(x, y, r) is then gov- 
erned by (see Appendix) 



OA 
97 



— = -(l + iA)A + £ + iT\A\ 2 A + i5V\A, (8) 



where is the transverse Laplace operator. The pa- 
rameter A in Eq. (8) is the detuning and is related 
to the linear phase accumulated by the light during 
one round-trip, ip = — 27rA/J r . The driving field is 
taken to be Ei n = 2ir£ / T exp(— iu>ot) and the coeffi- 
cient of the nonlinear term describes the combined ef- 
fect of the nonlincarities in both material layers: T = 
J 7 (tr Ir + 7iJl) /(27r), where the subscripts R and L de- 
note material parameters in the RHM and lhm, respec- 
tively. A self-focusing nonlinearity (r > 0) is assumed in 
the remainder of this paper; the case of a self-defocusing 
medium can be studied in a similar way. Finally, the 
diffraction coefficient S is given by 
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(9) 



Diffraction in the ring cavity depicted in Fig. 1 acts 
with the usual sign in the rhm, and can be compen- 
sated for by negative diffraction in the lhm. In order 
to get a proper understanding of the dynamics arising 
from the presence of the lhm layer in the cavity, we im- 
plement the procedure of Lugiato and Lefever to build 
a mean-field model (or LL-model) describing the evolu- 
tion of the field in the cavity [35]. This model is based 
on three additional assumptions about the optical cavity. 
First, we consider the RHM and the lhm to be adequately 
impedance matched, so that reflections at their interfaces 
can be neglected. As a result, the backward field (Ej,) is 
negligible. Note that this condition could also be im- 
posed in the system, by adding an isolator in the ring 
cavity. Second, we assume that the aspect ratio of the 
cavity is very large (high Fresnel number). In this case, 
the boundary conditions in both transverse directions do 
not affect the central part of the beam. Third, we con- 
sider cavities that are shorter than the diffraction and 
the nonlinearity space scales. 

For ease of notation, the time coordinate t has been 
normalized with respect to the lifetime of a photon in 
the cavity. This results in the slow dimensionless time 



Equation (8) is similar to the well-known LL-modcl. 
However, there is an essential difference: the diffraction 
coefficient can become negative. From Eq. (9), we deduce 
that this occurs when the length of the LHM layer is 
larger than —k-^ln/k-R,- We will demonstrate below that 
this negative diffraction may have dramatic consequences 
on the pattern formation process. 



IV. LINEAR STABILITY ANALYSIS 

The homogeneous steady state solutions (hss) A s of 
Eq. (8) are £ = [1 + i(A - T |A s | 2 )]yl s . The intracavity 
intensity \A S \ 2 as a function of |£| 2 is single- valued for 
A < y/3 and multiple- valued for A > \/2>. We performed 
the stability analysis of these HSS. With transverse pe- 
riodic boundaries, the deviation from the steady state is 
taken proportional to exp (i k • r + Ar) with k = (k x , k y ) 
and r = (x,y). Setting A = leads to a fc-dependent 
marginal stability condition: 



1+ (A + 5k 2 -2T\A s \ 2 y 



r 2 u s 



(10) 
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The condition dX/dk = then yields the expression of 
the threshold associated with the modulational instabil- 
ity (often called Turing instability): \A C \ = l/y/T and 
r|£ c | 2 = A 2 - 2 A + 2. At this bifurcation point, the 
critical wavelength is 



A c = 2tt 



' 7rJ"(A; L Z R + fc R j L ) 
k R k L (2 - A) 



(11) 



If the thickness of the lhm layer is small, Zl < 
—khln/kR, the system will have a positive diffraction co- 
efficient (6 > 0) and will behave like any traditional non- 
linear ring cavity. The more interesting situation occurs 
when Zl > — ^l^rAr- In that case, the diffraction coeffi- 
cient is negative, i.e., 5 < 0. The hss are unaltered, but 
their stability is completely different from the positive 
diffraction case. The upper branch is now stable what- 
ever the value of A. This means that the presence of a 
lhm completely inhibits the Ml in the monostable regime 
(0 < A < a/3). In contrast, for A > 2, an additional 
portion of the lower branch exhibits Ml, i.e., when 1/T < 



\A\ 2 < |A_| 2 , 
intensities |A 



or equivalently, when £ 

|2 



<£<£-. The 



(2A ± V A 2 - 3)/(3r) correspond ing 
to the driving fields £± = ^/\A±\ 2 [l + (A - |A±| 2 ) 2 ] are 
the coordinates on the bistability curve of the two limit 
points of the homogeneous hysteresis cycle. The behav- 
ior of the negative diffraction resonator is illustrated in 
Fig. 2 for A > 2. In Fig. 2(a), we have plotted the HSS 
as a function of the driving field; the domain of insta- 
bility is indicated by dashed lines. Figure 2(b) shows 
the marginal stability curve. In Fig. 2(c), the thresh- 
old associated with the MI as a function of the detuning 
parameter is plotted. From the latter figure, it is clear 
that the stability domain is bounded and that only a 
small portion of the lower HSS branch is affected by the 
mi. This is in strong contrast with traditional resonators 
where the upper HSS branch is unstable and where the 
instability domain is unbounded. 




FIG. 2: (a) Bistability curve of the intracavity intensity 
(| A\ 2 ) vs driving intensity (\£ | 2 ). (b) Marginal stability curve. 
The domain of instability is filled, (c) Intracavity inten- 
sity vs detuning (A). The plane is divided by the curves 
A = A ± and A = A c , with |A C | 2 = 1/r and |A±| 2 = 
(2A ± VA 2 - 3) /(3r). The domain of MI is filled. 
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V. DISCUSSION 

When the driving field amplitude is in the range £ c < 
£ < £- and A > 2, a finite band of linearly unstable 
Fourier modes k exist [see Fig. 2(b)], which trigger the 
spontaneous evolution of the intracavity field amplitude 
towards a stationary, spatially periodic distribution that 
occupies the whole cavity space available in the trans- 
verse directions. An example of a dissipative structure in 
the one-dimensional system obtained by numerical simu- 
lation of Eq. (8) is shown in Fig. 3. The initial condition 
consists of some small random noise added to the HSS 
and periodic boundary conditions are assumed. To com- 
pare with the analytical results, we calculate the wave- 
length (A c = 3.14) given by the linear stability analysis 
[Eq. (If)] and the wavelength of the stable dissipative 
structure obtained numerically (A num = 3.23). A very 
good agreement between the two wavelengths is obtained. 



FIG. 3: Stable dissipative structure emerging in the one- 
dimensional case when the resonator is biased near the thresh- 
old of mi. The real and imaginary parts of A are plotted with 
solid and dotted lines, respectively. Parameters are A = 10, 
6 = -2, T = 1 and £ = 9.25. 



A weakly nonlinear analysis performed in the vicinity 
of the threshold for mi shows that the branch of peri- 
odic solutions always appears subcritically for A > ||. 
Changing the diffraction coefficient in these simulations 
will simply change the scaling of the pattern, as can be 
seen from Eq. (8) by nondimensionalizing x and y us- 
ing the transformation (x,y) = V6(X,Y). Since we can 
make the diffraction coefficient of the resonator under 
study as small as we desire by choosing the thicknesses 
of the two material layers appropriately, we can in princi- 
ple produce dissipative structures with sizes smaller than 
the diffraction-limited size \flT~l (27rfc ). The combined 
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use of LHMs and rhms therefore enables us to circum- 
vent the diffraction limit, which must be compared to 
the sub-diffraction-limited imaging of lhm plates [36]. 




FIG. 4: Dynamics of the pattern formation in the two- 
dimensional case, (a)-(b) Initial pattern formation, (c)-(d) 
Hexagonal structures emerge, (e) Hexagonal structures are 
replaced by stripes, (f) Self-organized spots appear. The 
labels correspond to the points (a)-(f) on Fig. 5. The param- 
eters are the same as those in Fig. 3. Maxima are plain white 
and the integration mesh is 256 x 256. 

We carried out similar simulations in a square-shaped 
domain using periodic boundary conditions. The two- 
dimensional patterns are obtained by perturbing initially 
the HSS in the same parameter range as in Fig. 3. The 
following sequence of transverse patterns is obtained: 
hexagons, stripes, and self-organized spots. This behav- 
ior is shown in Fig. 4(a)-(f), where we have plotted the 
real part of the intracavity field during the spatiotcmpo- 
ral evolution. From Fig. 5, we see that the average inten- 
sity grows exponentially, which is clearly due to the am- 
plification of unstable wavelengths around A c = 2n/k c , 
as can be expected from the linear stability analysis. Fi- 
nally, the field amplitude reaches the high HSS. Note 
that the transition between the two HSS is abrupt, but 
takes place after a fairly long transient close the lower 
HSS branch. By varying the input field amplitude, we 
found that this phenomenon occurs for all input fields in 
the range £ c < £ < £_. This spontaneous spatiotempo- 
ral up-switching process therefore leads to the truncation 
of the homogeneous hysteresis cycle. 

To summarize, we presented a description of light prop- 
agation in a ring cavity filled with a left-handed and a 
right-handed material, both materials exhibiting a Kerr 
type nonlinearity. We derived a mean-field model for the 
spatiotemporal evolution of this resonator and showed 
that its diffraction coefficient can become negative if the 
thickness of the lhm layer is large enough. In this case, 
the LHM stabilizes the upper branch of the homogeneous 
steady state curve, and destabilizes a small part of the 
lower branch for high detuning (A > 2). In the latter pa- 
rameter range, we demonstrated that pattern formation 
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FIG. 5: Time evolution during the two-dimensional up- 
switching process described. The average maximum value of 
the real part of the intracavity field is plotted. 

can take place. However, a spontaneous two-dimensional 
up-switching process dominates the spatiotemporal dy- 
namics over the formation of stable patterns, resulting in 
the truncation of the homogeneous hysteresis cycle. 
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Appendix: Mean-field model 

The mean-field approximation has already been ap- 
plied to systems with boundary conditions [37]. Here, 
the propagation model is different, and we detail how it 
is possible to simplify the system of four coupled prop- 
agation equations, and the boundary conditions to only 
one equation of evolution, similar to the equation of the 
LL-model [35]. 

1. Hypothesis and notations 

We consider bidirectional propagation in a ring cav- 
ity with impedance matching at the interfaces between 
the two layers of classical and left-handed materials. In 
order to neglect the backward field with respect to the 
forward field in certain expressions, we also assume that 
the reflections at the other interfaces are limited. In cases 
where this assumption would not be verified, it is possi- 
ble to add a circulator in the cavity, so that the condition 
holds. 

Below, the letters R, L, and C will refer to the field, 
nonlinearity and dispersion in, respectively , the right- 
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handed medium, the left-handed medium, and the cavity 
outside of the two layers. The quantities A and B are 
(resp.) the envelope of the forward and the backward 
fields travelling in the ring cavity. 

With theses notations, for example, Ac(— Ir) is the 
electric field outside the right-handed medium, in the en- 
trance plane, while ^4r(— Ir) is the same field, at the same 
interface, but inside the medium. 

The mean-field approximation consists in approximat- 
ing the solution of the evolution in each layer, and for 
each forward or backward component as follows. In or- 
der to simplify the derivation, we rewrite Eq. (7) in the 
form 



3. Approximate equations of evolution 



8 A 

— = i[V+M{A,B)] A, 



(A.l) 



where I? is a linear differential operator, and Af(A, B) a 
nonlinear operator, i.e. 



V 



1 

2k 



t y2 _ ^ " 



2 dt 2: 



(A.2) 
(A.3) 



where a is the nonlinear coupling coefficient between the 
forward and the backward waves. 



2. Boundary conditions 

The geometry of the ring cavity leads to the following 
boundary conditions, 

AcHr) = t m E m +r m e lip A c (h), (A.4) 

Ar(-Ir) = t R Ac(-l R )+r R B R (-l R ), (A.5) 

A L (0) = A R (0), (A.6) 

A c (Il) = t L A L (l L )+r m r L e^ Bc(-Ir). (A.7) 

These expressions define the transmission and reflection 
coefficients at the interfaces. The phase factor e ttp ac- 
counts for the total phase accumulated in the ring cavity, 
outside of the two layers. 

Similar relations can be written for the backward field, 

BcHr) = t R B a (-l a )+r R Ac(-ln), (A.8) 

B L (0) = B R (0), (A.9) 

B l (Il) = t L B c (l L ) + r L A L (l L ), (A.IO) 

B c (Il) = r in e^B c (-l R ). (A.ll) 



Assuming that 5( 
Eq. (A.I) leads to 



0(1), a first order integration of 



A(( + 5() = A(0+i8([V+M(A,B)]\ ( A + O(2) 
= e iSC[D+Ar(A,B)]\ ( a + 0(2) 
= e iSCT>\ <e iSCM(A,B)\ ( a + o(2) 
= e i8(M(A,B)\ (e iSCD\ ( £ + 0(2). 

(A.12) 

Now, we make the assumption that the reflection coeffi- 
cients are small quantities (tl = 0(1) — tr), and that 
the backward field is initially very weak (B = 0(1)) in 
comparison with the forward field, so that 



AcHr) 


= t in E in +r in e lv A c (h), 


(A.I3) 


Ar(-Ir) 


- t R Ac(-l R )+0(2), 


(A.14) 


MO) 


= A R (0), 


(A.I5) 


A c (Il) 


= t L A L (l L ) + 0(2), 


(A.I6) 


BcHr) 


- O(0)Br(-I r ) + O(1)A c (- 


-i R )(k.n) 


Bl(0) 


= Br(0), 


(A.I8) 


B L (l L ) 


= O(0)B c (l L ) + O(l)A L (l L ) 


, (A.I9) 


B c (h) 


= 0(1)B c (-Ir). 


(A.20) 



In addition, the first order integration of the propagation 
equation Eq. (A.12) shows that 



Ar(0) 
A L (l L ) 
Br(-Ir) 
Bl(0) 



O(0)Ar(-I r ), 
O(0)A L (0), 
O(0)B R (0), 
O(0)B L (l L ). 



From Eqs. (A.13)-(A.20), we deduce that 

A = O(0)E in , 
B = 0(1)A. 



(A.21) 
(A.22) 
(A.23) 
(A.24) 
(A.25) 



(A.26) 
(A.27) 
(A.28) 



This shows that the backward field will stay a first order 
quantity if it is initially so with respect to the forward 
field, and that we can neglect the linear coupling with the 
backward field in the boundary conditions Eqs. (A.13)- 
(A.16), or Eqs. (A.4)-(A.7), as the coupling is of second 
order. 

The nonlinear coupling between the forward and the 
backward field is incoherent, meaning that it depends 
only on the intensity of the fields, and involves no ex- 
change of energy between the two counterpropagating 
envelopes. This allows to relate, at first order, the for- 
ward and backward intensities everywhere in the cavity 
to the forward input intensity. 
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A reduction of Eqs. (A.4)-(A.7) and (A.12) provides 
A {n) {-l R )= (A.29) 

tin E vn + T in t L t R e lip X 



e il L V L e il R T> Re il R N R A (u-1)(_ 1r ^ 



which expresses how the envelope of the forward field is 
related to both this same envelope one roundtrip before 
and the driving field E in . 

This discrete infinite set of equations can be reduced 
to a single equation, by introducing the roundtrip time 
T, and approximating the discrete variations by a partial 
derivative with respect to the slow time t — nT: 



A (n) 



A(n-1) ~ T ^A 

dr 



r—n T 



We find directly that 
OA 



T— = tin E in + (pe ir -1) A, 



(A.30) 



(A.31) 



where V = I l T>l + IrPr + 1 lNl + IrNr accounts for the 
total propagation, and p — r in t L t R e t{p . 

The backward field is proportional to the forward field, 
with \B c (-Ir)\ = \r R + r L t L t 2 R \\Ac{-l R )\ + 0(2), and 
similar expressions for \Br(—Ir)\ = Or\Ar(— Ir)\, and 
|-Bl(0)| = 0l|Ar(O)|. The exact expression of Or and 9l 
is not important to complete our derivation. Below, it is 
only needed to remember that 9 = 0(1), and that 9 is 
real. 

The linear relation between the forward and the back- 
ward intensities simplifies the calculation of the nonlinear 
operators, as 

AT(A, B) = 7 (\A\ 2 + a\B\ 2 ) = 7 (l + o9 2 ) \A\ 2 , (A.32) 

By redefining 7' = (l + o9 2 ) 7, we can completely hide 
the existence of the backward field, that appears as a sec- 
ond order correction to the nonlinear coefficient. There- 
fore 



II 
k L 



£)1 



2,L 



Ir&2,R. 



+ (Ili' l + Iri'r)\A\ 2 . 



dt 2 

(A.33) 



4. Lugiato-Lefever equation 



Together, equations (A.31) and (A.33) fully describe 
the dynamics of the cavity. In this section, we will assume 
that we are working near resonance of the cavity, so as to 
obtain the standard form of the LL-model, after proper 
normalization of the coefficients. 



Near resonance, p — \p\e 1 ^ is real, and ip accounts for 
the detuning with respect to the resonance (ip rts 0). In 
this particular case, 

pe iV = \p\e l ^ +v) = \p\ (1 + iip + iV) + 0(2). (A.34) 
Introducing this expression in Eq. (A.31), we find 

OA 

T— = t in E in + [(|p| - 1) + i\p\yj + i\p\V] A. (A.35) 

This equation is similar to the LL-model, which appears 
more clearly, if we redefine 



/ 1 - IpL 

t = — — — r, 



I" \P\ 



\P\ 



1 \P\ (II , I 



21-|p|V*L k R )' 
\P\ 



1-\P\ 
\P\ 



(l L /32,L + 1 R P2, R ) , 

(l L (l + a9 2 L ) lL + l R (l + a9 2 R ) 7 R), 



(A.36) 
(A.37) 

(A.38) 

(A.39) 

(A.40) 



(A.41) 



as in this case, the dynamics in the cavity near resonance 
is described by 



F)A B r) 2 A 

Z-=£ in -(l + zA)A + i5V\A + ^T\A\ 2 A. 

(A.42) 

When the system is driven by a continuous wave, the 
dispersive term disappears and we obtain Eq. (8). It is 
convenient to introduce the finesse of the cavity T = 
2tt jjqLj . If we neglect the second-order quantities 9 and 

assume that \p\ w \r in \ = t in = 1/V2, Eqs. (A.36)-(A.41) 
lead to the simple form 



, _ V2n 
T 

f = — E- 
A = -—ib, 

47T \k L k R 
T 



(A.43) 
(A.44) 
(A.45) 
(A.46) 



T = (IlIl + IrIr) ■ (A.47) 

Z7T 
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